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Abstract 

Identifying integrable coupled nonlinear ordinary differential equations (ODEs) of dissipative 
type and deducing their general solutions are some of the challenging tasks in nonlinear dynamics. 
In this paper we undertake these problems and unearth two classes of integrable coupled nonlinear 
ODEs of arbitrary order. To achieve these goals we introduce suitable nonlocal transformations 
in certain linear ODEs and generate the coupled nonlinear ODEs. In particular, we show that 
the problem of solving these classes of coupled nonlinear ODEs of any order, effectively reduces 
to solving a single first order nonlinear ODE. We then describe a procedure to derive explicit 
general solutions for the identified integrable coupled ODEs, when the above mentioned first order 
nonlinear ODE reduces to a Bernoulli equation. The equations which we generate and solve include 
the two coupled versions of modified Emden equations (in second order) , coupled versions of Chazy 
equations (in third order) and their variants, higher dimensional coupled Ricatti and Abel chains 
as well as a new integrable chain and higher order equations. 
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I. INTRODUCTION 



During the past few years considerable interest has been shown towards identifying a dif- 
ferential sequence of ordinary differential equations (ODEs) and studying their integrability 

and symmetry properties 1 ' 2 . An integrable differential sequence is a set of ODEs of which 

th 

each equation in the sequence is integrable. Here by an integrable n order ODE we mean 
that an equation for which the general solution x = f(t, C\, . . . , c n ), where q, i — 1, 2, . . . , n, 
are arbitrary constants, can be explicitly found (to within quadratures). 

To identify such a sequence of equations a general algorithm has been proposed by three 
of the present authors with Kundu 3 . The algorithm not only yields the well known Ricatti 
and Abel chains but also produces chains 4 ' 5 which involve inverse polynomials and more 
general systems. On the other hand, several recent studies have been confined to generation 
and analysis of only Ricatti and Abel chains 2 . 

The Ricatti chain is defined by D^z = 0, m — 0, 1, 2, 3, where = ^ + kx, k G 1R so 
that the action of 3 R generates the following family of differential equations 2 , namely 



and so on. The second member in this family is the Ricatti equation and the third member 
is the modified Emden equation and the fourth member is one of the subcases of the Chazy 
equation and so on. 

On the other hand considering the differential operator in the form B>a = ^ + kx 2 , fcet, 
one can generate a hierarchy of higher order Abel equations 2 , namely 



X 
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x = x + /cx 3 



x = x + Akx 2 x + & 2 x 5 , 



x = x + 5kx 2 x + 8kxx 2 + 9k 2 x 4 x + k 3 x 7 , 



(2) 



etc. The second member is a special case of the Bernoulli equation and the third member is a 
generalized van-der Pol oscillator equation and the fourth member is a subcase of the Chazy 
equation. It has been shown that all the equations in each of the hierarchy posses certain 
common properties that make them interesting both from physical and mathematical points 
of view. For example, as far as the Ricatti chain is concerned (i) all the equations in this 
chain admit a maximal number of Lie point symmetries, (ii) every equation of order m in 
this hierarchy can be linearized and transformed into a linear ODE of order m + 1 and (hi) 
the dimensional reduction of a linear equation of order m + 1 leads to the Ricatti equation 
of order m. As far as the Abel chain is concerned all the equations in this hierarchy posses 
a family of first integrals J t k that depend on time as a polynomial. In addition to the above, 
specific equations in both the categories have also been investigated in detail. In particular, 
the second equation in the Ricatti chain, namely the modified Emden equation, has been 
a central attraction for mathematicians and physicists for more than a century (for more 
details see Refs.[6-10]). 

As we mentioned earlier, in Ref. [3] it has been demonstrated that one can generate 
both the Ricatti and Abel chains of equations from the damped linear harmonic oscillator 
equation, U + C\U + c^U = 0, where c\ and c 2 are arbitrary parameters, and higher order 
linear ODEs by introducing a nonlocal transformation, U = x a e^^ xm+ ' y ^dt, where a and 
m are constants and j3 and 7 are arbitrary functions of t. Substituting the above nonlocal 
transformation in the damped harmonic oscillator equation and restricting the parameters 
appropriately (a — 1, (3 = k, 7 = 0, C\, c 2 = 0, m = 1) one can get the third member in the 
Ricatti chain. Similarly, taking a = 1, f3 = k, 7 = 0, ci, c 2 = 0, m = 2 one can get the third 
member in the Abel chain. Now considering a third order linear ODE U+ciU+czU +C3U = 0, 
where c±, c 2 , c 3 are constants, and substituting the above nonlocal transformation in this 
equation, and restricting the parameters suitably, one can get the fourth member of both 
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the chains. The authors have also presented a straightforward method of finding explicit 
general solution for all the equations belonging to these chains. 

Even though several works have been dedicated to explore the physical and mathematical 
properties of these two chains no attempt has been made to identify/generate a sequence of 
coupled integrable ODEs and construct their general solutions, which are in general difficult 
tasks. The aim of this paper is to make some progress in these directions and present 
some of the important results regarding how to generate two different classes of integrable 
coupled ODEs and the method of constructing their general solutions. We show that the 
problem of solving the chain of coupled integrable ODEs generated by this procedure can 
be ultimately reduced to solving a single first order nonlinear ODE. In other words, the 
problem of constructing the solution of the system of nonlinear coupled ODEs effectively 
reduce to deducing the solution of the equivalent first order nonlinear ODE. We deduce the 
solution of this first order equation for the parametric choice for which it reduces to an 
integrable Bernoulli equation. We find that for specific choices of system parameters the 
system of coupled nonlinear ODEs exhibit isochronous oscillations, which is an interesting 
area of study in its own merit 11 . In the process of identifying two new integrable class of 
equations we have obtained a higher dimensional generalization of the Ricatti and Abel 
chains. The two dimensional Ricatti chain has been identified as 




(3) 
(4) 
(5) 



xi = 0, 
x 2 = 0. 

x\ + (a\Xi + a 2 x 2 )x\ = 0, 
x 2 + {aiXi + a 2 x 2 )x 2 = 0. 

X\ + 2(aia;i + a 2 x 2 )x\ + [a\X\ + a^^x^ + (aiXi + a2X2) 2 xi = 0, 
x 2 + 2(diXi + a 2 x 2 )x 2 + (aiXi + a 2 x 2 )x 2 + (01X1 + a 2 x 2 ) 2 x 2 = 0. 

x'i + (3(aiXi + a 2 x 2 ))xi + x(aiXi + a 2 x 2 ) + 3(aiii + a 2 x 2 )xi 

+ (aiXi + a 2 x 2 )(3xi(aiXi + a 2 x 2 ) + 3xi(aiXi + a 2 x 2 ) + (aiXx + a 2 x 2 ) 2 Xi) = 0. 
x 2 + (3(oiXi + a 2 x 2 ))x 2 + 2/(01% + a 2 x 2 ) + 3(oix a + a 2 x 2 )x2 

+ (aiXi + a 2 X2)(3x 2 (aiXi + a 2 x 2 ) + 3x 2 (aiXi + a 2 x 2 ) + (a^i + a 2 x 2 ) 2 x 2 ) = 0, 

and so on, where W R — (4r + (a±xi + a 2 x 2 )Y . This two dimensional Ricatti chain can be 
further generalized to arbitrary dimensions and the details are given in Sec. I VI I The two 
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dimensional Abel chain is given as 

K \ " I Xl = 0, (7) 

x 2 = 0. 



(8) 
(9) 



X\ + (aixf + i^^D^i = 0' 
x 2 + (aixf + a 2 x\)X2 = 0. 

Xi + (2(a,\x\ + g^x^Ai + X\((a\x\ + a 2 x 2 ) 2 + 2aiXiiri + 202^2X2) = 0, 
X2 + (2(aixf + a 2 x 2 ))x 2 + x 2 ((aixf + a 2 x\) 2 + 2a\X\X\ + 202X2X2) = 0. 

x'i + 3(aixf + a 2 x\)'±\ + 6(aiXiXi + a 2 x 2 x 2 )(xi + x 1 (aixf + ^2^2)) 
+2xi(aiXi + a 2 x 2 x 2 ) + 2aiX\x\ + 3(aixf + a 2 x 2 ) 2 + x 1 {d\x\ + a 2 xl) 3 = 0, 
x 2 + 3(aix 2 + a 2 x 2 )x 2 + 6(diXiXi + a 2 X2X 2 )(x2 + x 2 (aix 2 + a 2 x|)) 
2x 2 (ciiXi + a 2 x 2 x 2 ) + 2a2X2x| + 3(aix 2 + a2x|) 2 + x 2 (a\x\ + a 2 xl) 3 = 0> 

and so on, where = (jj + (a\x\ + 02X3))'. Some of the interesting equations in the above 
chains are the coupled modified Emden equations (j3j), coupled generalization of Chazy type 
equation and the coupled generalized Duffing- van der Pol oscillator equations ([9]), the 
details of which are given in Sees. II III and |VJ The n-dimensional versions of the above two 
chains are in Sec. IVI1 

In addition to the n-dimensional Ricatti and Abel chains, we have also identified a new 
integrable chain. The simplified two dimensional version of this new integrable chain is given 

as 





a\x\x 2 = 0, 
x 2 + a\X\x\ = 0. 

+ 2aiXix 2 xi + a\x\x\ + ai(xiX2 + xii 2 ) = 0, 
'x 2 + 2aiXiX2X 2 + a\x\x\ + ai(xiX2 + X1X2) = 0. 

x'x + 2ai(xiX2Xi + xiXix 2 + x 2 X2) + ai(3xfx 2 xi + x\x 2 ) 
+ai(xix 2 + Xix 2 + 2iix 2 ) = 0, 

x 2 + 2ai(xix 2 x 2 + x 2 Xix 2 + xix 2 ,) + a\ (3x1X3X2 + Xix|) 
+ai(xix 2 + xix 2 + 2iix 2 ) = 0, 



(11) 
(12) 

(13) 



(14) 
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and so on, where Ojy = (A + a^XiX^) . The above integrable chain can be further generalized 
to n-dimensions and the details are also given in Sec. I VI I 

The plan of the paper is as follows. In Sec. [Ill we point out the nonlocal connection 
between two uncoupled damped linear harmonic oscillators and certain coupled nonlinear 
oscillator equations. We also present a method of identifying suitable nonlocal transforma- 
tions for this purpose. We then discuss the method of constructing the general solution for 
the identified set of nonlinear two coupled ODEs in Sec. IIHI Generalizing this analysis, 
in Sec. \IV\ we consider a system of n uncoupled damped harmonic oscillator equations 
and generate a class of n coupled nonlinear second order ODEs. In Sec. |V] we extend the 
procedure to identify a class of integrable third order coupled nonlinear ODEs. To make our 
studies systematic here also we first consider a set of two coupled third order linear equations 
and by introducing the same nonlocal transformations we identify a class two coupled third 
order nonlinear ODEs. This procedure is then extended to n-coupled third order ODEs. 
In Sec. I VI I we show that the results obtained for second order and third order coupled 
nonlinear ODEs can be extended to arbitrary I order coupled nonlinear ODEs. As special 
cases, in addition to the n-dimensional generalization of Ricatti and Abel chains, we identify 
a new n- dimensional integrable chain of nonlinear ODEs. In Sec. IVHI we consider a system 
of uncoupled inhomogeneous linear equations with variable coefficients and briefly discuss 
the applicability of our procedure. We summarize our results in Sec. Villi 

II. NONLOCAL CONNECTION BETWEEN SECOND ORDER LINEAR ODES 
AND COUPLED NONLINEAR ODES 

While studing the dynamics of a set of two coupled modified Emden equations 12 , 

x + 2(fax + k 2 y)x + (fax + k 2 y)x + (fax + k 2 y) 2 x + X±x = 0, 

y + 2(kiX + k 2 y)y + (fax + k 2 y)y + (fax + k 2 y) 2 y + X 2 y = 0, (15) 

we have identified the fact that the system (JT5l) can be transformed into two uncoupled 
harmonic oscillator equations, namely U + X±U = and V + X 2 V = 0, where Ai,A2 are 
constants, through the nonlocal transformation U = xe^ klX+k2V ^ dt and V = ye^ klX+k2V ^ dt . 
Interestingly the solution for Eq. (TL5]) can be derived from the known solution of the un- 
coupled linear equations U + XiU = and V + X 2 V = 0, say U = a(t) and V = b(t). This 
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can be done by noting the fact that x and y are also related with U and V through the set 
of first order coupled nonlinear ODEs of the form 



U V 
x = —x - (fax + k 2 y)x, y = —y - (fax + k 2 y)y. 



(16) 



Since jj = ^ and y = ~ are known functions, Eqs. (JTBjl can be integrated straightforwardly 
and rearranging the constants suitably one can deduce the general solution for equation ( fl5l) 
in the form 



x(t) 

y(t) 



Asm(uit + Si) 



1 - ^ cos^t + 8 1 )-££ cos(u 2 t + 5 2 

B sin(u 2 t + 5 2 ) 
1 - 5 cosM + - ^ cos(a; 2 t + <5 2 ) ' 



Ai, w 2 



+ 



tu 2 



A 2 
< 1 



(17) 



where A, 5\, 5 2 are four arbitrary constants. 



A. Generalization 

Now the question naturally arises as to whether one can generalize the nonlocal transfor- 
mations and identify a chain of coupled nonlinear ODEs as in the case of the scalar case 3,13 . 
A systematic investigation reveals the fact that this can indeed be done. 

For this purpose, we consider a set of two uncoupled damped linear harmonic oscilla- 
tors, defined by the system of uncoupled second order linear ordinary differential equations 
(ODEs) 

U + c n U + c l2 U = 0, V + c 21 V + c 22 V = 0, (18) 

where Cn, Ci 2 , c 2 \ and c 22 are arbitrary parameters for the present. However, they can be as 
well functions of t as we point out later in Sec. IVIIII Introducing a nonlocal transformation, 

U = x a e Jf(x,y,t)dt^ y = y f> e fg(x,y,t)dt^ ^gj 



where / and g are two arbitrary functions of their arguments, in ([18]) we obtain a set of two 
coupled second order nonlinear ODEs of the form 

. 2 

x + (a-l)— + (2/ + cn)i + - \f + c n f + c 12 + f]=Q, (20a) 
x a l J 

. 2 

y + (0 - l)^- + (2g + c 21 )y + f [g 2 + c 21 g + c 22 + g] = 0, (20b) 

y p 



; df df .df .df , . dg dg . dg . dg _ , . , 
where / = — = — + x— + y— and g = — = — + x— + y—. Obviously choosing 
at at ox ay dt at ox oy 

a = (3 = 1 and the functions / and g to be linear in x and y, Eq. (120]) reduces to (|T5|) . 

In order to find a family of integrable nonlinear ODEs in the class (120]) which are con- 
nected with the damped linear harmonic oscillator equations ([18]) through the nonlocal 
transformation ( fl9]) and deduce the general solution for the coupled nonlinear ODEs from 
the solutions of the linear ODEs (1181) . we find that the task essentially consists of two parts. 
First one has to choose specific forms of / and g and fix the target coupled nonlinear ODEs 
and then one can look for a suitable method to derive the general solution of the nonlinear 
ODEs. In the following we describe a procedure which successfully takes care of both the 
steps. 

To start with we observe that the nonlocal transformation (TT9T) yields a system of two 
first order coupled nonlinear, nonautonomous ODEs of the form 



x 

x = — 
a 



jj - f(x,y,t) 



y 



y -g(x,y,t) 



(21) 



where U and V are the solutions of the linear ODEs (Tl8]) . One may note that since jj and y 
are some known functions of l f the first terms on the right hand sides of both the equations 
have explicit time dependent factors. 

Note that for the choice f(x,y,t) = an(t)x + ai 2 (t)y and g(x,y,t) = a 2 i(t)x + a 2 2(t)y, 
Eq. ( ]2T]) reduces to the time dependent two dimensional Lotka-Volterra (LV) equation 

x = x(a n (t)x + a 12 (t)y + &i(t)), y = y(a 21 (t)x + a 22 {t)y + b 2 {t)), (22) 

where bi(i) = jj and b 2 (t) = y. The parameters a, (3 have been absorbed into bi and a^-, 
i,j = 1,2. Studies have been carried out on the dynamics of this system 14 ^ 16 . However, 
we find that Eq. ( 122]) is not integrable for arbitrary choice of the coefficients a^-, and bi, 
i,j = 1,2, but for special forms they are integrable as we see below. 

In general it is easier to integrate Eq. ( 12 ip than to solve Eq. ( 12 Op and obtain x(t) and 
y(t). The question now boils down to the fact for what forms of / and g Eq. f[2~T]) can be 
integrated to yield the general solution. While analysing the form of Eq. ( |2T]) we find that 
upon choosing / and g in a symmetric form a first integral can be obtained by suitably 
rewriting and integrating the resultant equations. Of course the form of / and g have to be 
fixed in such a way that the resultant set of coupled first order ODEs falls into a coupled 
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Bernoulli family of equations so that the general solution can be obtained. To be specific, 
we identify two specific forms of / and g. They are 

(i) / = ( ai (t)x p + a 2 (t)yi) m + h(t), g = ( ai (t)x p + a 2 (t)y^) m + b 2 (t) (23) 

with the parametric restriction — = q (details of this parametric restriction are given in 
Sec. IHip and 

N N 

(n) f=(52ai(t)x"y«) m + h(t), g = ($>(t)x«j,«) m + (24) 
i=i i=i 

with the parametric restrictions — + = k, i = 1,2, . . . , N. Here a^s and b^s, i = 
1, 2, . . . , N, j — 1, 2, are arbitrary functions of t, and p, q, pi, qi are real numbers satisfying 
the above parametric restrictions and A; is a constant. For both the cases we find that the 
resultant Bernoulli equation can be integrated to yield the solution which in turn leads to the 
general solution of the nonlinear coupled second order ODEs with appropriate redefinition 
of the integration constants as we see in the following sections. 



III. INTEGRABLE TWO COUPLED SECOND ORDER NONLINEAR ODES 



In this section we illustrate the procedure for two coupled second order ODEs. The two 
different forms of / and g given by Eqs. ( l23l) and (124]) are treated separately as Case 1 and 
Case 2 in the following. 
Case : 1 

For the forms of / and g given in ( )23|) . Eq. (I2T]) reduces to a system of coupled Bernoulli 
type equations, 



(25a) 
(25b) 



X 


u 


a 


u 


y 


~v 


p 


V 



Substituting this form of / and g into Eq. ( 1201) we obtain the corresponding set of coupled 
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nonlinear second order ODEs as 



x + (a - 1)^- + [2((a lX p + a 2 y q ) m + b x ) + c u ] x + - 



((a lX p + a 2 y q ) m + h) 2 + c 12 
+c u {(a 1 x p + a 2 y q ) m + b 1 )+m {a x x p + a 2 y q ) m ~ l ^px^x + a 2 qy q - l ij 



+a lt x p + a 2t y c - 



nt 



(26a) 



x 2 



y+(P-l)^- + [2((a lX p + a 2 y q ) m + b 2 ) + c 21 ] y + V - 



(( ai x p + a 2 y0) m + b 2 ) 2 + c 22 



+c 21 (( ai x p + a 2 y q ) m + b 2 )+m (a lX p + a^)" 1 ' 1 { a.px^x + a^y^y 



+a lt x p + a 2t y q } + b 



>2t 



0. 



(26b) 



Solution of Eq. (1261) can be written down by solving the coupled Bernoulli type equations 
( 1251) provided the parametric condition ^ = q is satisfied . For this purpose, we multiply 
the first equation of (T2"5]) by ay and the second by fix and subtract the latter from the former 
to obtain 

'u V 



ayx — (3xy 



+ b 2 (t)-h(t) xy. 



(27) 



U V 

Note that the general form of the solutions of the linear ODEs ( fl8l) can be written down as 

U(t) = he mit + I 2 e m2t , V(t) = he mzt + he mi \ (28) 

where I±, I 2 , I3, and I± are four arbitrary constants, mi and m 2 are the solutions of the 
auxiliary equation m 2 +cnm+ci 2 = 0, while m 3 and are the solutions of m 2 +c 2 im+c 22 = 
0. Now dividing both sides of ( 1271) by xy and integrating the resultant equation ( 12T1) one 

gets 



x(t) 



v y 



(29) 



We also observe that the integration constant has been absorbed with the arbitrary constants 



in ^. Substituting ( 1291) in the second equation of Eq. (1251) we get 



V \ y y 



ai 



V 



y°<+ a 2 y 



The above equation reduces to the Bernoulli type equation 

y = f(t)y - g(t)y mq+1 



(30) 



(31) 
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for the choice — — q. Note that for the simple choice a — fl, we have the relation p = q. 
Integrating now equation (}30l) we obtain y(t) and once y(t) is known, x(t) can be obtained 
straightforwardly through (129|) (where again the integration constant has been absorbed into 
that of U and V), 



x(t) 



y(t) 



i -i- m i 

1 "T /3 



eft 



aqm 



1 1 r/6 2 rft 



-i i mq 



Oi 



(32a) 



(32b) 



Note that the general solution contains four arbitrary constants through U and V, see Eq. 

(EHD. 

As a specific example of (126]) . let us consider for the simple choice a = (3 = p = q = m = l 
and ax, a 2 , bx, b 2 are constants. For this parametric choice Eq. (126]) reduces to 

x + 2(axx + aiy + bx)x + (axx + a2y)x + (axx + a^y + bx) 2 x 

+cxx(x + x(axx + a 2 y + h)) + c 12 x = 0, 

y + 2(axx + a 2 y + b 2 )y + (axx + a 2 y)y + (axx + a 2 y + b 2 ) 2 y 

+c 2 i(y + y(axx + a 2 y + b 2 )) + c 22 y = 0. (33) 



The general solution of Eq. (!33j) reads 



x(t) 



Ue 



-bit 



y(t) 



1 + / (axUe-^ + a 2 Ve- b *)dt' 
Ve~ b2t 



(34) 



1 + J(axUe- b ^ + a 2 Ve~ b2t )dt' 

As one expects the general solution of the nonlinear system depends upon the solution 
of the damped linear harmonic oscillator equation too. As a consequence the solution of 
Eq. fl33|) may be periodic or decaying or growing type depending upon the nature of the 
parameters, Cn, c 12 , c 2l , and c 22 which appear in (ITS]) . We note that for the parametric 
choice bx, b 2 , Cn, c 2 x = 0, Eq. (1331 becomes Eq. f[T5"j) and admits isochronous oscillations 
where the frequency is independent of amplitude or initial conditions 11,12 , see Eq. (|T71) . 
Also note that for the choice Cn = c 12 = c 2 i = c 22 = and bx = b 2 = 0, Eq. (|33l) reduces to 
the third member of the two coupled Ricatti chain (jSJ). 
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Case : 2 



Apart from the above discussed form of / and g which reduces Eq. ([21]) to an in- 
tegrable system of coupled Bernoulli type equations, we find the following form, / = 

f N \ m , N sm 

^OityaPy*] +b 1 {t)&ndg= l^OityaFy*] + b 2 (t), where Pi) q h i = 1, 2, . . . , N, 

\i=l J \i=l J 

and m are arbitrary real constants, also reduces Eq. (12 ip to the following coupled Bernoulli 
type equations 



x 



x 



a 



U 
U 



JV 



,i=i 



P 



v 

V 



N 



h(t)- [J2<t> Pi y 9i 



,i=i 



.(35) 



However we find that the general solution of Eq. (j35l) can be deduced for the specific 
parametric restrictions ^ + <&= constant, where i = 1, 2, . . . , N. Substituting this form of 
/ and g in Eq. ( 120]) we get the associated integrable coupled nonlinear second order ODEs, 



x+(a-l)— + 2x[(%2 aiX Pi y q T + h + c u 



N 



i=l 



X 

a 



JV 



i=l 



+ C12 



JV 



N 



1=1 



i=i 



JV 



+a it x 1H y 



nt 



y](aiPiX Pi l y qi x + aiqiX Px y q * x y 



(36a) 



* i=l 



JV 



JV 



JV 



+c 2 i (X^VT + &2 +m(J2 a i xPi y q T~ 1 



8=1 



i=l 



+ 6: 



2/ 



0. 



J^a^y*) +6 2 +c 22 

.i=l / 
r V 

^(a^x p '~ V l i + a i q i x n y qi '' 1 y 

(36b) 



i=i 



Again, multiplying the first equation of (J3"5]) by ay and the second of ([35]) by /fa and 
subtracting the latter from the former one obtains an equation for (ayx — f3xy) in exactly 
the same form as that of Eq. (I27p . As before, dividing the resultant equation on both sides 
by xy and integrating one gets x(t) = ^e^ b2 ~ bl ^ dt y^ a . Substituting this in the second 
equation of Eq. ( 135]) we get 

JV 



,i=i 



V 



(37) 
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The above equation reduces to the Bernoulli type equation ( 13T]) for the following parametric 
restrictions 



Pifi i il ^ • 1 O AT 

a m 



(38) 



where k is a constant. The parametric restrictions f )38|) contains equations as i runs from 

TV 

1 to N and this set of conditions implies that all the terms in the expression aiX Pi y 9z are 

i=i 

of the same degree. The parametric restrictions ( 138]) reduce equation (|37|) to 







r 
v 



A 



i=l 



f(b 2 -bi)dt 



(39) 



Integrating equation ( |37j) we obtain y(£) and once y is known, x can be deduced straight- 
forwardly through (1291) . Doing so one obtains 

U^e-^f bldt 



x(t) 



1 + ^ / F^e^" 



(fc-i) (i-fc) 



A' 



/&2<it 



v(*) 



i=l 



C7 



,/(62-6i)rft 



Pi\ ™ 



-,(40a) 



i(fe-i) 



1 1 fb 2 dt 



1 + ^ / V^e^r 



(fc-i) (i-fc) 



N 

E 

.i=i 



/(6 2 -6i)dt 



Pi \ m 



-. (40b) 



As an interesting example, for the choice a = /3 = m = pi = qi = N = l, Eq. ( l36i) reduces 
to 



x + 2x{a\xy + b\ + en) + a\x 2 y + alx 3 y 2 + x 2 y(a u + 2a\b\ 

+x(b 1 c n + b u + C12 + b?) = 0, 

y + 2y(a l xy + b 2 + c 21 ) + a x x 2 y + a 2 x 3 ?/ 2 + x 2 y(a u + 2aib 2 ) 

+x(b 2 c 21 + b 2t + c 22 + bj) = 0. 

The general solution of the above equation is 

U e -f b i dt 



(41a) 



(41b) 



x(t) 



y(t) 



[1 + 2 J ciiUVe-f^+^dt] 1 

y e -fb 2 dt 

[l + 2j ai UVe-f^+ b ^ dt dt]' 



(42a) 
(42b) 
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Note that solution (1421) can become isochronous for the choice b\ = b 2 = and suitable 
choice of other parameters in f l4"TT) . Also note that for the choice en = c\ 2 = C21 = C22 = 0, 
b\ = b 2 = and ai, a 2 are constants, Eq. ( l4"Tj) reduces to the third member of the new 
integrable chain f lT3"j) . 



IV. INTEGRABLE ARBITRARY n-COUPLED SECOND ORDER NONLINEAR 
ODES 

The above analysis can be generalized in principle without much difficulty to a system of 
arbitrary n-coupled second order ODEs and classes of solvable ones from the linear ODEs 
can be identified. In order to do so we first investigated the case of three-coupled second 
order ODEs. Then from the results of two- and three-coupled second order ODEs, we have 
generalized the results to n-coupled second order ODEs. For brevity, we are not presenting 
here the results of the three-coupled second order ODEs (one can refer to Ref.[17] for details), 
but straightaway present the results for a system of n-coupled ODEs. 

Let us consider a system of n uncoupled linear second order ODEs of the form 

Ui + CnUi + c i2 Ui = 0, i = l,...,n, (43) 

whose general solution can be expressed as Ui(t) = Iue mut + l2ie m2it , where In and I 2 i, 
i = 1, 2, . . . , n, are arbitrary constants and m^'s and m^s are the solutions of the auxiliary 
equations. 
Case : 1 



Now introducing the nonlocal transformations analogous to ( |23i) 

Ut = ^e/^^r^WK, i = l,...,n, (44) 

in f T45|) , one can transform the uncoupled damped linear harmonic oscillators into the sytem 
of n-coupled nonlinear second order ODEs of the form 

. 2 n r n / n \ 

+ ( ai - 1)^ + 2({j2 y + + — (J2 ^r -1 ( w* -1 ^ + ^ ) 

X i , a i L i=1 \ i=1 / 



3 

r / 1 1 



+ - 



n \ 2 n 

a 3 xf) m + bA + Cil ((^ a jX p /) m + k) + 

i=l ' 1=1 



c t2 



+ c a x = 0, % = 1,2, ...n. (45) 
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Now we confine our attention on obtaining the general solution for Eq. ( f4"5]) which can be 
done by following the same methodology described earlier for the two coupled case. Then 
Eqs. (JHJ can be brought to the form of a set of coupled Bernoulli type first order ODEs, 

. Ui 



- ((^2 a jX f) m + bi)xi, i = l,...,n. 



(46) 



One can reduce the above system of 72,-coupled first order ODEs to a single first order 
equation as was done in the case of n = 2 (vide Eq. [301) . For a better understanding we 
demonstrate the procedure for n = 3. For this case the above system of equations ( )46|) 
reduces to 



Xi 

X\ = — 



x 2 

%2 = 

a 2 



x 3 

x 3 = — 

a 3 



^ - 61 - (a^ 1 + a 2 xf + a 3 xf) m 



El _ b2 _ ( aiX W + a2X P2 + a3X P^m 
U2 



u 3 

— 63 - (aixf + a 2 x p 2 2 + a 3 Xg 

(-4 



P3 \m 



(47a) 
(47b) 
(47c) 



We solve the above system of equations using the same method used to obtain ( |29l) . Mul- 
tiplying (I47a[) by aix 3 and (I47cl) by ot 3 x\ and subtracting the latter from the former one 
obtains 



ot\X 3 Xi — a 3 xix 3 



JJ- ~ JJ- + h(t) - h(t) j .r s .r->. 



(48) 



Dividing both sides by X\X 3 and integrating the resultant equation fj48l . after absorbing the 
constant of integration, one gets 



El P J(h-bx)dt T a 3 

u 3 e 



(49) 



Multiplying (147b|) by a 2 x 3 and (I47c[) by 03X2 and subtracting the latter from the former one 
obtains 



U 2 U 3 \ 
"2 ■'■■.',■>■■■> - <^i- 2 .r : ; = I — - — + b 3 (t) - b 2 (t) I X 2 X 3 . 

Dividing both sides by x 2 x 3 and integrating the resultant equation (1501) one gets 



(50) 



x 2 



El e .f(b3-b2)dt x a s ^ " J 
^3 



(51) 
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Substituting (gSJ) and ([EE]) in Eq. (H7cj) we get 



2:3 

X3 = — 

a 3 



lh 



U 



pi 

a, Pl a 3 



lh 



a 2 P2 a 3 

a 2 1 ~JP3 



(52) 



In a similar way the system of n coupled Bernoulli type equations (|46j) can be brought to 
the single first order equation for the nth variable x n (t), 



X ■;■ 



a,, 



Un 



— b„ — 



i=l 



TT \ ~ p j°' n 



(53) 



Utilizing the same idea given in the previous case and using the parametric restriction 
^aT = ^c^T = '" = one can inductively deduce the general solution of Eq. (jlSjl as 



Xi(t) 



-J- fbidt 

1 (2 OLA " L 



■y 1 mp„ 

On 



7M* 



_, a 7l 1 

mp n a i 



dt 



i = l, 2,..., to. (54) 



The general solution contains 2n arbitrary constants, coming from that of solutions of Eq. 
(PI). 

Case : 2 

Substituting the nonlocal transformation U = £^ e ^ ( ^= iaj(i)n ^ lX ^ )m+felW , in Eq. (JS} we 
get the set of n coupled second order nonlinear ODEs, 



N n 

+ («. _ 1)^ + 2X, ( (£ a, J] ) m + fe, ) + 

i=i fc=i 



TV 



V n 

i=i fe=i 



x 



N n 

E%n 

j=l k=X 



X 



Pjk\m—X 



+ 



i=l I fc=i 

£«iII*?T+& 



(=1 



\ 



/ 



3=1 k=X 



H (c,2 + &it, 



a. 



Af n 

E°dl 

3=1 k=X 



X 



Pjk\m 



Cixx = 0. 



(55) 



Proceeding as before, for the parametric choice 
solution of the above equation is given by 



i=X 



^ + 1 = -, j = 1,2, ... ,N, the 
a,- m 



Xi(t) 



l £> OLA J 1 



1 + ^1 s{t){U n e-I bndt )^dt 

. Here the set Ui, i = 1, 2, 3, . . . , n, is the 



an 1 



(56) 



where s(t) 



EN T—rn— 1 

i=1 ai[[ j=1 



u 



4 ^siK-b 3 )dt 

e 3 



solution of the linear system of ODEs ( I43p . 
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V. INTEGRABLE COUPLED THIRD ORDER NONLINEAR ODES 



A. Two coupled third order ODEs 

Having studied the nonlocal connection between second order linear and nonlinear ODEs 
now we generalize the results to a system of third order ODEs. For this purpose let us 
consider the following system of uncoupled third order linear ODEs, 

U + c u U + c 12 U + c 13 U = 0, V + c 2 iV + c 22 V + c 23 V = 0, (57) 

where % = 1,2, j = 1,2, 3, are arbitrary parameters, and the same nonlocal transforma- 
tion ffT9]) . 

JJ _ x <x e J f(x,y,t)dt y _ y/3 e f g(x,y,t)dt ^ 

Substituting the above nonlocal transformation in fl57|) we obtain the following system of 
two coupled third order nonlinear ODEs, 

x + (3(a - 1)- + 3/ + c n + -f x )x + -f y y + (a - l)(a - 2)^ + (^—^(3/ + c n ) + 3/, 

■ -LrU- + (3f y + -f xy )xy + (3/ 2 + 2c n / + c 12 + -./',, + — //. + + 3f t )x 

a a a a a 

+(-fyt + —ff v + —f y )y + -(Sfft + f + c u f + c n ft + fu + c 13 + c 12 f) = 0, (58a) 
a a a a 



y + (3(/3 -l) V - + 3g + C2i + \g y )y + ^g x x + (p - l)(/3 - 2)^ + (^— ^(3<? + cai) + 3g y 
V pa y 2 y 



+ p9yy)y + {3g x + -j9xy)xy + {3g + 2c 2 ig + c 22 + — g yt + —gg y + —9y + 3gt)V 
+(-jj-f*t + ~^99x + ^j-g x )x + ~(3ggt + g 3 + c n g 2 + c 21 g t + g a + c 23 + c 22 g) = 0. (58b) 

As in the case of the system of second order ODEs, we find that there are two sets of / and 
g (vide Eqs. (|23|) and (|2^|) ). for which the above equation becomes integrable. We briefly 
discuss below the general solution of Eqs. (158]) for both the choices. 
Case : 1 

Substituting / = { ai {t)x p + a 2 (t)y q ) m + 6i(t) and g = (ai(t)x p + a 2 {t)y q ) m + b 2 {t) in (pD we 
get a system of two coupled third order nonlinear ODEs. The solution of this set of nonlinear 
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ODEs can again be obtained by solving the following coupled Bernoulli type equations 

U 



x 



x 
a 

y_ 





-b 2 (t)-( ai (t)x p + a 2 (t)yT 



(59a) 
(59b) 



which is obtained through the nonlocal transformation U = x a e^^ ai ^ xF+a2 ^ yq ^ m+bl ^ dt , 

y = y a e J(( ai (t)xP+a 2 (t)y"r+b 2 (t))dt j n (jggj) thg aboye forms of f and g 

We find that the above equation (15 9p is exactly the same as that of equation (125]) whose 
solution has been discussed in Sec IHIj except for the difference that U and V are now 
solutions of the third order linear ODEs ( 15 7p . Therefore we directly write down the solution 
of equation (159]) (vide Eq. ( 132]) ) as 



x(t) 



y(t) 



1 I rag 



U] { Ue-f bldt y +a 2 [Ve- J' b2dt 



3.1 m 



dt 



amq 



1 -±fb 2 dt 



1 + 7 



J ai (Ue-fWy + a 2 (ve-f'^ dt y 



dt 



1 1 
mq 



(60a) 



(60b) 



where U and V are now the general solutions of the system of third order linear ODEs (jo 
which can be given in the form 



U(t) = he mit + I 2 e m2t + / 3 e m3 \ V(t) = he mit + J 5 e m5 * + he 



Jfn&t 



(61) 



where 7i, I 2 , 13, ■ ■ ■ , Iq, are six arbitrary constants and mi,m 2 , are the solutions of the 
auxiliary equation 



m 3 + Cum 2 + c\ 2 m + c 13 = 0, 



while m 4 , m 5 , m 6 are the solutions of the auxiliary equation 



m 3 + c 2 im 2 + c 22 m + c 23 = 0. 



(62a) 



(62b) 



We note that for the parametric choice Cn = C13 = c 2 i = c 2 % = the general solution 
of Eq. (1571) . U and V, become periodic and for this parametric choice with b\ = b 2 = 0, 
we find that the solution (l60j) exhibits isochronous behaviour. As specific example let us 
consider the choice a — /3 = p — q — 1, b\ — b 2 — 0, ai, a 2 constants. Then the nonlocal 
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transformation U = x <* e S(*i(t> p +<v(t)y q ) m +b 1 (t)dt ^ y _ y a e J(a 1 (t)xP+a 2 (t)yT+b 2 (t)dt re d nces to 
the form U = X e^ aix+a2V ^ dt and V = y e ^ aiX+a2y ^ dt . Substituting this form of nonlocal 
transformation in 

U = 0, V = (63) 
we get the following set of third order nonlinear ODEs, 

x + 3(aii + a 2 y)x + 2(a x x + a 2 y)x + (a\x + a 2 y)x + 2(a x x + a 2 y)x + (a x x + a 2 y) 2 x 
+ (aix + a 2 y){x + 2{a x x + a 2 y)x + (a\X + a 2 y)x + {a\x + a 2 y) 2 x) = 0, 

V + 3(aix + a 2 y)y + 2{a x x + a 2 y)y + {a x x + a 2 y)y + 2(aix + a 2 y)y + {a x x + a 2 y) 2 y 
+(a 1 x + a 2 y)(y + 2{a x x + a 2 y)y + {a x x + a 2 y)y + (a x x + a 2 y) 2 y) = 0. (64) 

Equation f[6"l"|) is nothing but a system of two coupled Chazy equation XII (with N = 2 and 
parametric restrictions A = 0, B = in Ref. [18]), which has been studied in detail in Refs. 
[18-23]. 

/ N \ m / N \ m 

Case : 2 We substitute / = ( ai(t)x Pi y qi J + bx(t) and g = ( ai(t)x Pi y qi I + 6a(*) 
(which is the same as used earlier in Sec. II III for second order ODEs) to obtain 



i 3 



x + 8i(x,y,t)x + S 2 (x,y,t)y + (a - l)(a - 2)— + 8 3 (x,y,t)x + 5 4 (x,y,t)y 



+5 5 (x, y, t)xy + 8 6 (x, y, t)x + 5 7 (x, y, t)y + 5 8 (x, y, t) = 0, (65a) 
y + e x (x, y, t)x + e 2 (x, y, t)y + ((3 - l)((3 - 2)— + e 3 (x, y, t)x 2 + e 4 (x, y, t)y 2 

y 

+e 5 (x, y, t)xy + e 6 (x, y, t)x + e 7 (x, y, t)y + e 8 (x, y, t) = 0, (65b) 

where 5i{x, y, t) and €i(x, y, t), i = 1, 2, 8, are functions of the indicated variables and we 
do not explicitly present them here as they are complicated expressions and one may refer 
to Ref. [17] for details. 

As we have pointed out in the previous case, the solution of the two-coupled second order 
nonlinear equation (|36|) and the solution of the two-coupled third order nonlinear equation 
( 16"5"|) are the same except for the form of the solutions of the linear ODEs. Therefore one can 
write down the solution of Eq. ( 1651) also directly using the solution of the two-coupled second 
order nonlinear equation ( 136]) with the parametric restriction — +^+1 = — , i — 1, 2, . . . , N, 
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as 



x(t) 



U^e-^f bldt 



N 



1 • V / \—<^ lh -" l! [Y. a -(^ (j ' 



y(t) 



riLJl V j b. 2 dt 



. i=l 
1 1 



u 



Pi\ rn 



dt 



i 

»(fe-i) 



-,(66a) 



fc-i 



V I s e f> 



Jb 2 dt 



N 

E 

i=l 



L7 

», I — e 



/(6 2 -6i)d* 



— , (66b) 



where U and are now the general solutions of the system of third order linear ODEs ([5 
given as in Eq. ( l62l ). 



B. Arbitrary n-coupled third order ODEs 

We conclude this section by extending the ideas to a sytem n-coupled third order nonlinear 
ODEs by extending the ideas of the previous discussion. Let us consider a system of n 
uncoupled linear third order ODEs of the form 



Ui + c a Ui + c i2 Ui + c i3 Ui = 0, i = 1,2, ...,n. 



(67) 



Now introducing the nonlocal transformation 



one can derive the following set of coupled nonlinear ODEs, 



x\ + ( 3(<*i - 1)- + 3/ + ax j i % + V ™ + 



j'=i 



3/ j 3/ 

3 n 



fa,- - 1} 



,1 , 



(3/i + Cii) + 3/ iXi x 



•2 



^_ ^ f X^ J ^ ^ ^2iX ' y 



* t=i 



3=1 J 



+(3tf + i,.f, + c i2 + 3/«)±, + V^/,, + —J)/,.,- , + ^)x 3 + ^(3fJ it + ff 

' ™- a,- a,- a,- 



(68) 



+Cji/? + Cji/ft + fat + Cffl/< + Cjs) = 0, i = 1, 2, . . . , n. 



(69) 



We again consider two different forms of fiS for which one is able to deduce the solution of 
Eq. ( 169]) separately. 
Case : 1 
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Considering the form of fi as fi = (VJ™ cijX^ 3 ) m + bi) and imposing the parametric restric- 
tions = 2^ = ... — p n , we obtain the following solution, as in the case of the coupled 
second order ODEs given by Eq. fl54l) . 



Xi(t) 



tt ola — — f bidt 



1 + 



nip,, 



3=1 



-fbjdt) «n 



dt 



SZ- , i = l,2,...,n, (70) 



where U^s are the solutions of the system of third order linear ODEs (IBTl) containing 3iV 
arbitrary constants. 
Case : 2 



Considering the form of fi as fi = {Y^f=i a j{t) n«=i x K K ) m + h{t) and imposing the 



parametric restrictions ^ u - + 1 = — , j = 1, 2, . . . , N, we obtain the following solution 

^— f a,- m 



TTCtj — — f bidt 



1 + 



A ~ 7 s(t)(U N e-f bndt )^dt 



-, i = l,2,...,ra, (71) 



where s(£) 



i=l a UI 7 =i 



4 e ^ /(&»-&*)*" 



and CZj's are the solutions of the system 



of third order linear ODEs (JS7J). Again (ITT]) is of the same form as (156]) obtained for the 
case of coupled second order ODEs 



VI. INTEGRABLE COUPLED / tn ORDER NONLINEAR ODES 

Having studied the nonlocal connection that exists between linear and nonlinear ODEs 

for the second and third orders we finally build a theory which is applicable for arbitrary 
th 

I order equations. 

Consider a system of two uncoupled linear ODEs 

{ d l d (l -V d \ 

(tf + *^ + -" + <*-'5 + < *) lr M = ' (72) 
where q/s, i = 1,2, j = 1,2,... I, are arbitrary constants. The nonlocal transformation 
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(|T9|) connects ( 172]) to the set of coupled nonlinear ODEs of the form 

Df ] + c n Dt 1] + ■■■ + cu-iD^ + di)x = 0, 



\D® + c 21 Dt 1} + ■■■ + cn-iD™ + c 2 ^j y = 0, (73) 

where = (ai a_1 (| + xf(x,y,t))) and D% = (/3?/ /3 ~ 1 (^ + yg(x, y, t))) 1 . As we have 
seen, the solution of Eq. (!73|) can be deduced from the nonlocal transformation and the 
solution of the system of linear ODEs ( 1721) only for specific forms of f(x, y, t) and g(x, y, t), 
which are given below. 
Case : 1 



Solution for Eq. ( 173]) can be deduced by the procedure developed earlier for the choice 
f{x,y,t) = (ai(t)x p + a 2 y q ) m + h(t) and g(x,y,t) = (ai(t)x p + a 2 y q ) m + b 2 (t) with the 
parametric restriction ^ = q. The solution for this case is the same as the one given by Eq. 
(132]) . except for the fact U and V are now solutions of the fi 1 order linear ODEs ( 1721) which 
can be given as 

i i 
U(t) = J * e "^ V(t) = ( 74 ) 



i=l i=l 



where U and Ii, i = 1,2, ... ,1 are 21 arbitrary constants while rrii and rhi are respectively 
the / roots of the auxiliary equations 

m! + cum'" 1 + . . . + c n = and rh! + c 2 im^ 1 + . . . + c 2i = 0. (75) 

Case : 2 

For the choice f(x,y,t) = (y%=i a^x^Y + b 1 {t) and g(x, y, t) = (j^ti <k{t)x^y^ m + 
b 2 (t) and with the parametric restriction + qi + 1 — —, i — 1, 2, . . . , N, Eq. (I73|) can be 
solved. The solution for this case is also of the same form as given by Eq. (HOl) . except that 
now U and V are solutions of the Eqs. (172]) containing 2/ arbitrary constants as given in 

flzU). 

th 

The results can be straightforwardly extended to the n-coupled I order ODEs. We do 

th 

not discuss here separately the three coupled I order ODEs. 
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A. Arbitrary n-coupled / order ODEs 



Consider a set of n uncoupled linear /th order ODEs of the form 
/ d l dP -1 ) d \ 

U +c "^ + -" +c -* +c «j c, '" ) = - (76) 

where Ci/s, i — 1, 2, . . . n, j — 1, 2, . . . I, are arbitrary constants. The nonlocal transformation 
Ui = x^eJ ' fi( x i> x 2>-> x n> t ) connects ( 1761) to the nonlinear ODE of the form 

(dP + caDf-V + ... + c^D? + cuj x l = 0, i = 1, 2, . . . , n, (77) 

where Df ] = (o^ -1 ^ + x i f i (x 1 ,x 2 , . . . , x n , t))) 1 . 

Again we find that the solution of Eq. (176"]) can be obtained for the following two cases. 
Case : 1 

fi = (Ylk a k x k k ) m + h{t)) with the parametric restriction = = ... = p n leads us 

to the solution in a form exactly similar to Eq. ( 15^1) for the arbitrary n-coupled third order 

th 

ODEs except that U^s i — 1, 2, . . . , n, are now solutions of the uncoupled system of I order 
linear ODEs ([76]) as 

i 

U t (t) = J2^ m ^ (78) 

3=1 

where 1^, i = 1, 2, . . . , n and j = 1,2,...,/ are arbitrary constants and solutions of 

the auxiliary equations 

m^ + Ciim^ 1 + ... + Cjj = 0, « = 1,2, . . . , n. (79) 

We note here that for the special choice of parameters en = = . . . = cu = 0, oti = 1 and 
Pk = m = 1, k — 1,2, ... ,n, Eq. ( ]77j) reduces to the n- dimensional coupled Ricatti chain 
whose explicit form can be given as 

d n V 

17 + ^ Qfc^fc j Xi = 0, i = 1,2, . ..,7i. (80) 
fe=i / 

For the choice n = 2 the above equation reduces to the two dimensional Ricatti chain given 
inEqs. 

Similarly for the parametric choice cu — — . . . = cu — 0, p/. = 2, k — 1, 2, . . . , n, and 
m = ttj = 1, Eq. (I77|) reduces to the coupled n-dimensional Abel chain whose explicit form 
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can be given as 



n x ' 



^ + ^2a k x 2 k \ x { = 0, i = l,2,...,ra. (81) 

k=l / 

Note that the above equation reduces to the two dimensional Abel chain given in Eqs. 
(T7])-(TTU]) for the choice n = 2. 

n 

Case : 2 fa = (X)^Li a jW rifc=i xV k h ) m + with the parametric restrictions, _|_ 

i=i 

1 = — , j = 1,2, . . . , N, gives us again the solutions xAt) in exactly the same form as 
m 



given in Eq. (J56J) for the arbitrary n-coupled third order equations except now the set Ui, 
i = 1, 2, 3, . . . , n in (1561) is the solution of the uncoupled linear system of ODEs (!76|) as given 
in (JTSD. 

In addition to the above n-dimensional Ricatti and Abel chains (vide (l80l) and ( JHTT) ) a 
new ro- dimensional integrable chain can be identified which can obtained from Eq. (1771) with 
the parametric choice en — c& = . . . — c# = 0, ctj = 1 and it reads 

N 

+ r' 



^ + £ a ^Il^T| ^ = °- t = l,2,...,n, (82) 



i=i k=i 

and J2i=iPij + 1 = ^, J = 1,2, ...,iV. By choosing N — 1, pjk = 1 and n = 1, one can 
obtain Eqs. (jTT])- (fl4]) . 



VII. INHOMOGENEOUS LINEAR ODES AND ASSOCIATED COUPLED NON- 
LINEAR ODES 

In our studies, on connection between uncoupled linear ODEs and coupled nonlinear 
ODEs through nonlocal transformations, we have assumed all along that the set of uncoupled 
linear ODEs are homogeneous with constant coefficients. However, this is not mandatory 
for our procedure to hold good. For example, we can choose a system of I order uncoupled 
linear ODEs of the form 

d l U- d l ^ 1 U- dU- 

~df + Cli ^~dF^ + " ' + Q ~ M ^"df + Cl ^ Ui = ' = 2 ' ' ' ' ' n ' (83) 

where Cjt(t), j — 1, 2, . . . I, i — 1, 2, . . . , n and Fi(t) are functions of It is well known that 
the general solution of ( 18H|) can be written as 

Ui (t) =a u (t) U u (t) + a 2l (t) U 2i (*) + ...+ a K (t) U h (t) , (84) 
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where Z7ji(i)'s, j = 1, 2, ... Z for a given z (= 1, 2, ... n) are Z linearly independent solutions 
of the homogeneous part of Eq. fl83l) and the coefficient functions ajj(t)'s can be given by 
quadratures which involve I arbitrary constants, for a fixed i. So as long as the general 
solution of (|83|) is known, then the entire procedure developed in Sees. HTllVIl completely 
goes through. The corresponding coupled nonlinear ODEs will have coefficients involving 
Cji(t)'s which will be now functions of t. 

VIII. FINAL COMMENTS 

In this paper, we have developed a novel method of identifying two classes of integrable 
coupled nonlinear ODEs of any order from linear uncoupled ODEs of the same order by 
introducing suitable nonlocal transformations in the latter. We found that the problem of 
solving these classes of coupled nonlinear ODEs of any order, effectively reduces to solving a 
single first order nonlinear ODE and we have deduced the general solution for the parametric 
choice for which this first order nonlinear ODE reduces to Bernoulli equation. For suitable 
choice of parameters we find that the coupled nonlinear ODEs can also exhibit isochronous 
behaviour. 

In this paper we have focussed our attention only on identifying coupled integrable equa- 
tions through specific type of nonlocal transformations connecting dependent variables only. 
There are several possible generalizations : (i) We have considered only two specific in- 
tegrable choices of Eq. (I2TI) . One can look for other possible integrable cases too. (ii) 
The functions /, g in the transformations can be functions of x, y also, (iii) The inde- 
pendent variable can also be transformed and (iv) One may even consider transformations 
connecting uncoupled integrable nonlinear ODEs and coupled nonlinear ODEs. All these 
possibilities can lead to several new classes of integrable coupled nonlinear ODEs. Some of 
these problems are being pursued at present. 
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